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Nonlinear dynamics of the interface of dielectric liquids in a 
strong electric field: Reduced equations of motion 
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The evolution of the interface between two ideal dielectric liquids in a strong 
vertical electric field is studied. It is found that a particular flow regime, for which 
the velocity potential and the electric field potential are linearly dependent functions, 
is possible if the ratio of the permittivities of liquids is inversely proportional to the 



ratio of their densities. The corresponding reduced equations for interface motion 
are derived. In the limit of small density ratio, these equations coincide with the 
well-known equations describing the Laplacian growth. 
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It is well known that the flat interface of two dielectric liquids is unstable in a sufficiently 
strongvertical electric field. The dispersion relation for the surface waves has the following 
form Q : 

Pi+p 2 47r(pi + p 2 ){£i + £2) P1 + P2 
where k is the wave number, u is the frequency, g is the acceleration of gravity, a is the 
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surface tension coefficient, p\ and p2 are the mass densities of lower and of upper liquids 
(Pi > P2), £\ and e 2 are the dielectric constants of fluids. The external electric field strengths 
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under and above the interface, E\ and E 2 , are related by the expression 

e l E l = e 2 E 2 . (1) 



It is seen from the dispersion relation that, if the electric field is sufficiently strong, 

£1 +£ 2 



EiE 2 > -— 2 —^/ga{pi-p 2 ), 

{Si — e 2 ) 
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the second term in right-hand side of the dispersion relation dominates for the waves with 
wave numbers in the range 

g(£i + e 2 )(pi ~p 2 ) k E l E 2 (e 1 -e 2 ) 2 
E 1 E 2 (e 1 -e 2 ) 2 a(ei + e a ) 

Then u 2 oc k 2 and, hence, we can separate the dispersion relation into two branches 



(±) _ ±ick c 2 = / 2 \ 

47r(p 1 +p 2 )(e 1 + e 2 )' 

For one branch, small periodic perturbations of the surface increase exponentially with the 
characteristic times (c/c)" 1 , while, for the other branch, these perturbations attenuate. In 
such a situation, we can restrict our consideration to the increasing branch = +ick, 
that essentially simplifies the problem of describing the evolution of the interface at the 
linear stage of the development of instability. The buildup of perturbations of the surface 
inevitably transforms the system to a state in which its evolution is determined by nonlinear 
processes. Then, in the general case, splitting into the branches becomes impossible. 

In this paper we will show that, for the particular case E\p\ = e 2 p 2 , we can extract the 
separate branches from the equations of motion. This makes it possible to reduce by half 
the number of equations required for describing the evolution of the boundary. The reduced 
equations coincide with the well-known equations describing the Laplacian growth in the 
limit of small ratio of liquid densities. An important point is that the Laplacian growth 
equations not only define a subclass of particular solutions of the problem, but they also 
describe the asymptotic behavior of the system. 

It should be noted that the behavior of the interface of two fluids in normal electric or 
magnetic field (these problems are similar from the mathematical point of view) is usually 
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investigated in the quasi-monochromatic approximation (see |2l m, li| and the references 
therein). This approach allows one to obtain immediately an equation for the complex 
amplitude of surface waves. However, the applicability of such an equation is limited by 
the condition of the smallness of the slopes of the surface. The development of instability 
can violate this condition. In the strong- field limit, the approach developed in the present 
work provides a way of studying the interface behavior at essentially nonlinear stages of 
instability development. 

Consider the evolution of the interface of two ideal liquids of infinite depth in an external 
vertical electric field. In the unperturbed state, the boundary of the liquid is a flat horizontal 



E\E 2 \£\ — £ 2 
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surface. Let the z axis of the Cartesian coordinate system is normal to the unperturbed 
interface. The function r)(x,y,t) specifies the shape of the deformed boundary, i.e., the 
liquids occupy the regions z < r)(x,y,t) and z > r)(x,y,t), respectively. It is convenient 
for the subsequent analysis to choose an origin of coordinates so that the level of liquids is 
determined by the expression z = —vt. In other words, the origin moves with respect to the 
interface at a certain constant velocity v. 

Let us assume that the motion of both liquids is potential. The velocity potentials for 
incompressible liquids $x and $ 2 satisfy the Laplace equations, 



V 2 $i = 0, 



0. 



(3) 



with the following conditions at the boundary and at infinity: 
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z = ri{x,y,t) 



z — ► — oo 
z — > +oo 



(V^i-V^ 2 ), z = r]{x,y,t), (4) 

(5) 

(6) 
(7) 



where y>i and <^2 are the electric-field potentials in and above the liquid, and d/dn denotes 
the derivative along the normal to the interface. The expression on the right-hand side of 
the dynamic boundary condition (nonstationary Bernoulli equation) is responsible for the 
electrostatic pressure at the interface between two ideal dielectric liquids in the absence 
of free electric charges ^7]. The evolution of the interface is determined by the kinematic 
relation, 



drj <9$i 



(Vl77-Vl$i), z = r}{x,y,t). 



at dz 

The electric potentials ipi and ip 2 satisfy the Laplace equations, 



o, vV 2 = o. 



(8) 



(9) 



Since the electric field potential and normal component of the displacement vector have to 
be continuous at the interface, we should add the following conditions at the boundary: 



£l 



dn 



¥>2, 



£2 



d(p 2 



dn ' 



z = r](x,y,t), 
z = r){x,y,t). 



(10) 
(11) 
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The system of equations is closed by the condition of the electric field uniformity at an 
infinite distance from the surface: 



^2 



-Eiz, 
-E 2 z, 



-oo, 



z — > +00. 



(12) 
(13) 



Let us show that a flow regime, wherein the harmonic potentials of velocity and of electric 
field are linearly dependent functions, is possible for certain relations between the problem 
parameters. Suppose that 



.1/2 



if 2 = &$ 2 (4vrp 2 /£2 



,1/2 



(14) 



where a and b are unknown constants. It is necessary to verify that the initial equations 
of motion (J3|)- fl~3|) are compatible with these relations. Substituting them into (J3J) and 
(P2JHHH!), we obtain 



Pi 



d$i (V$i) 2 

dt 2 



■P2 



<9$ 2 (V$ 2 

2" 



:a6(ei— e 2 ) 



$ 1 a(47rp 1 / £l ) 1 / 2 = $ 2 6(47rp 2 /e 2 ) 1 / 2 , 

b( P 2S2) 1/2 , 



a{pi£i) 1 



dn 



d<& 2 
dn 



y/Plp2 (VgrVgg) 
\/£i£2 2 

z = rj(x,y,t), 
z = rj(x,y,t), 



^ 1 ^-zE 1 a- 1 (A7T Pl /e 1 )- 1 / 2 , 
$ 2 -> -2E 2 r 1 (47rp 2 /£ 2 )- 1/2 , 



z 
z 



-00, 
-00. 



z = r)(x,y,t), 

(15) 
(16) 
(17) 

(18) 
(19) 



For the system of equations (JHJ), ©-© and (|THj) - (|Tn)l to be compatible (it is overdetermined 
in the general case), the conditions (J5)l— f7)l must coincide with the conditions ()17| ) -([19 p . and 
the condition (jHJ) must coincide with the condition ()15j) . 
It is apparent that the conditions © and (JT7J) coincide if 



a(p 1 e ] 



1 1/2 



1/2 



(20) 



In view of Eqs. ((TJ) and (|2*0jl. the conditions at infinity (JBJ), ((7J) and (fTSj). ([T9*]) are consistent 
if the auxiliary parameter v takes the following value: 



v — a 1 v , 



E 1 (4rrp 1 /e 1 )- 1/2 > 0. 
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Finally, we consider the condition under which the dynamic (|15|) and kinematic (JHJ) rela- 
tions coincide. Let us write Eq. (jHj) in the form which does not contain function rj explicitly. 
With the help of the formula (ffilj) . the boundary condition (fi7)j) can be rewritten as follows: 

e^&x = E 2 - l & 2 , z = r)(x,y,t). 

Differentiating this expression with respect to time or spatial variables, we arrive at 



drj 
dt ' 



dz dz 

d$i d§ 2 

dz dz 



z=r] 



Z=JJ 



~di dt 



z=rj 



-[eVJi- V ± $ 2 ] 2= 



where e = £2/^1 is the ratio of the permittivities. These relations allow us to eliminate 77 
from Eq. (j^J. We obtain from the kinematic boundary condition: 

<9$i <9$ 2 



dt 



dt 



-£(V$!) 2 + (V$2- V$i), z = ri(x,y,t) 



(21) 



Decomposing the velocities of fluids into the normal (<9$/<9n) and tangential (<9$/<9r) com- 
ponents in Eqs. (fT5|l and (J5TJ), and taking into account Eqs. (0), (fTfij) and pPjl. we get 



9$i <9$ 2 



[1 — p — a 2 e 1 +a 



dt 



dt 



~i 2 



<9f 



<9f 



dn 
1-e- 1 



[l — ps 2 — a 2 +a 2 e) 



1 2 



<9n 



~i 2 



0. 



<9r 

« = v( x ,y,t) 



0, z = r}(x,y,t) 



where p = pij P\- Clearly, these expressions can coincide only if the following conditions 
hold: 



p = e , 1 - pe - a z + a e = 0, 
1 - p - a 2 e- x + a 2 = 2 - 2e _1 . 

From this it is inferred that the equations are compatible provided that a 2 = 1 and also 



£iPi = s 2 p2- 



(22) 



The equation for the parameter a has two roots, = ±1, corresponding to different 

branches of solutions. 

Thus, we have proved that the functional relation (|14|) can be compatible with the equa- 
tions of motion if the condition (j22j) is valid. The corresponding flow regime is described by 
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the following equations: 

V 2 $x = 0, V 2 $ 2 = 0, (23) 

§ = ^ y/l + (V^) 2 , z = v(x,y,t), (24) 

( + \ fort 
^ = ^, z = T,{x,y,t), (25) 

p 1 Qi = p 2 &2, z = r)(x,y,t), (26) 
$i -> -a (= % z, z -> -oo, (27) 
$ 2 -> -a (= % <2, z -> +oo. (28) 

The reduction of the initial equations dSJ (H3J) to Eqs. (|2^j) - pH|) significantly simplifies the 
analysis of the interface motion. As will be discussed below, in the formal limit p-ij P\ — * 0, 
these equations describe the so-called Laplacian growth. 

Let us find the dispersion relation for Eqs. (J2HJ)-(J2BJ)- We will seek a solution in the form 



v = c ^-^t) _ a (±) V()tj 

where ci, C2, and Ci are small constants. These expressions correspond to a small-amplitude 
sinusoidal deformation of the initially plane liquid-liquid interface. After simple transfor- 
mations, we obtain the following relation between the frequency u and the wave number 
k: 

u, = iaW v k = fl W 1 = P2 \ El k. 

Pi + P2 V^Pl \Pl + PV 

It can be seen that, for the branch a = a < - + - ) = +1, initial perturbation will increase and, 
for a = = —1, it will attenuate. It should be noted that, with regard to Eqs. (0) and 
(I22|) . this expression coincides with the expression (J2J) specifying different branches of the 
dispersion relation for the unreduced equations of motion. 

Thus, if the condition (|2*2*|) is satisfied, the separation of two branches corresponding to 
solutions increasing and decreasing with time is possible not only in the linearized equations, 
but also in the initial nonlinear equations l{3* |) -()13| ) . 

The question arises as to whether the flow regime under consideration is stable. In 
other words, whether or not Eqs. (J25)l - (J25Jl describe the large-time asymptotic behavior 
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of the system. Stability of the increasing branch is evident at linear stages of the interface 
evolution, when the linearized equations of motion can be split into two independent systems. 
At the nonlinear stages, the equations do not split completely, and the stability problem 
becomes nontrivial. 

It turns out that the stability can be proved in the limiting case p\ ^> P2 and e\ <C e<i 
(the condition ()22|) can be violated). Then the evolution of the interface will be governed by 
the influence of the lower liquid. The equations determining the interface motion become: 

v 2 $! = o, vVi = o, 

ot 2 87rpi 

^ + (V$i • V</?i) = 0, z = n(x, y, t), 

(pi = 0, z = r](x,y,t), 

$! — > —VZ, Z — > — OO, 

ifi — > —EiZ, z — > — oo, 

where the kinematic boundary condition is given in the implicit form. If we introduce a pair 
of auxiliary potentials, 



these equations can be rewritten in the following symmetric form (compare with Refs. 

V 2^(±) = 0> ( 29 ) 

^+(V^±)) 2 = 0, z = V (x,y,t), (30) 

^r(+) = ^r(-) j Z = T){x,y,t), (31) 

¥ +) -> -v z, z -> -oo, (32) 
-»• z -> -oo. (33) 

Here, we set a = +1 and, as a consequence, v = v . 

One can readily see that these equations are compatible with the condition = 0, 
which corresponds to the situation of interest, where the velocity potential and the electric 
field potential are functionally related. For = 0, the set of equations (|2~9|l - (j53|l reduces 



8 



to 

V 2^(+) = 0> ( 34 ) 

On d^ {+) , 

¥ = ~d^~ ^ + (V±r/)2 ' * = 1&V>t)> (35) 

*W = 0, z = r]{x,y,t), (36) 

#(+) _ UoZ; ^ _ oc . (37) 

The same equations can be immediately obtained from Eqs. ()23j) - (j2"%|) in the limit p — > 0. 
They coincide with the equations describing the so-called Laplacian growth, viz., the motion 
of the phase boundary with a velocity directly proportional to the normal derivative of a 
certain harmonic scalar field (\E^ + ) in our case). Depending on the chosen frame of reference, 
this field may have the meaning of temperature (Stefan's problem in the quasi-stationary 
limit), electrostatic potential (electrolytic deposition), or pressure (flow through a porous 
medium). It is important for us that there are many known exact solutions to Eqs. f|34[)— 



(137)) . They describe the evolution of the inter 
dimples, and so on (see, for example, 



Ho, 



11 



ace up to the formation of "fingers" , cuspidal 
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Let us prove that the class of solutions of the motion equations (|2*5|) - (|3^|l corresponding 
to the reduced Eqs. (|51)l - (|S7j) is stable to small perturbations of potential It should be 

noted that the motion of the liquid-liquid boundary described by Eqs. (p^)) -(p?7 | is always 
directed towards the lower liquid; this is associated with the extremum principle for harmonic 
functions. Let function rj at the initial instant t be a single-valued function of variables x 
and y. In this case, for t > t , the inequality 

r](x, y, t) < T](x, y, t ) (38) 

holds for any x and y. This inequality remains valid for small perturbations of also, 
when the effect of potential \I/ <_ ' ) in relation (J3*T|) can be disregarded as compared to the effect 
of potential *ff( + \ and the motion of the boundary is described by the same Eqs. (|34 p -()37j ) . 

As regards the evolution of potential it is sufficient, for small |V\E ,( - _ ' I |, to consider 

the boundary condition ()30)) in the linear approximation. It takes the trivial form: 

*J _) = 0, z = r)(x,y,t). (39) 

This means that the potential does not change with time in the chosen reference frame (the 
origin moves relative to the interface with speed Vq). In the simplest case of a periodic 
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perturbation, the solution to Eqs. (j2HJ), and (|59*)l is given by 

vpH = Ae KZ sxn(nx), 

where k is the perturbation wave number, and A is a constant small amplitude. Let us 
denote the potential at the boundary z = rj by ip. We have 

ip(x, y, t) = (_) | z= „ = Ae«rt x ' v > t) sin (/car). 

Taking into account the inequality (J3HJ), we finally get 

\ij>(x,y,t)\ < \^(x,y,t )\ 

for any x and y at t > to, that is the value of the potential at the interface does 
not increase with time. Furthermore, since the level of the interface (the value of function 77 
averaged over the spatial variables) moves downwards at a constant velocity, it is evident that 
the potential \I>( _ ) relaxes to zero at the boundary. Thus, we have proved that Eqs. (|53jl - (|37jl 
describe the asymptotic behavior of the liquid-liquid interface in a strong vertical electric 
field. 

It should be noted that the results of this work can be used to describe the motion of 
the interface of two dielectric liquids in an applied electric field for other geometries of the 
problem. All one has to do is to modify the conditions (pj). (J7|), (|12j). and (JT3J). This will 
allow us to consider the interface dynamics in an oblique or tangential electric field, and also 
the dynamics of closed interfaces in an external field. 

In addition, the results of the above investigation can be extended to the case of two mag- 
netic fluids in a vertical magnetic field. For this purpose one should replace the electric fields 
E lf2 by the magnetic fields and the permittivities by the magnetic permeabilities 
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